Numerical study of spin quantum Hall transitions in superconductors with broken 

time-reversal symmetry 
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We present results of numerical studies of spin quantum Hall transitions in disordered supercon- 
ductors, in which the pairing order parameter breaks time-reversal symmetry. We focus mainly on 
p-wave superconductors in which one of the spin components is conserved. The transport properties 
of the system are studied by numerically diagonalizing pairing Hamiltonians on a lattice, and by cal- 
culating the Chern and Thouless numbers of the quasiparticle states. We find that in the presence of 
disorder, (spin-) current carrying states exist only at discrete critical energies in the thermodynamic 
limit, and the spin-quantum Hall transition driven by an external Zeeman field has the same critical 
behavior as the usual integer quantum Hall transition of non-interacting electrons. These critical 
energies merge and disappear as disorder strength increases, in a manner similar to those in lattice 
models for integer quantum Hall transition. 

PACS numbers: 74.40.+k, 73.43.Nq, 72.15.Rn 



I. INTRODUCTION 

Transport properties of quasiparticles in unconven- 
tional superconductors have been of strong interest to 
condensed matter physicists since the discovery of high 
T c cuprate superconductors. The cuprates, which are d- 
wave superconductors, support gapless nodal quasiparti- 
cle excitations; these nodal quasiparticles dominate heat 
and spin transport at low temperatures. Some heavy 
fcrmion superconductors are also known to have d-wave 
pairing. Another class of unconventional superconduc- 
tors that have been receiving increasing attention are 
p-wave superconductors. While the p-wave pairing was 
originally found in superfluid 3 Heji interest in it has been 
renewed more recently due to advances in two distinct 
systems. Firstly, fractional quantum Hall effects at fill- 
ings factor v = 5/22i are believed to be the manifesta- 
tion of the p-wave pairing of composite fermions^* 4 ***^ 
Secondly, a growing number of results on the unconven- 
tional superconductivity of Sr2Ru04 emerge in favor of 
a triplet-pairing order parameter^ 

In addition to the experimental relevance, unconven- 
tional, disordered superconductors are also of great in- 
terest for theoretical reasons, as they represent new sym- 
metry classes in disordered non-interacting fermion prob- 
lems that are not realized in metals. A classification of 
these symmetry classes have been advanced recentlyAifi 
Depending on the existence (or the lack) of time-reversal 
and spin-rotation symmetries, dirty superconductors can 
be classified into four symmetry classes, CI, Dili, C, and 
D in Cartan's classification scheme. These classes are be- 
lieved to complete the possible universality classes^ in 
disordered single-particle systemsAiS 

Systems of class C with broken time-reversal invari- 
ance but preserved spin-rotation symmetry exhibit uni- 
versal critical behavior which has been under extensive 
investigation recently>i2ii&iii5ii& Such systems, which 



can be realized in two-dimensional superconductors with 
d x 2_ y 2 + id xy symmetry, have a distinct signature of a 
critical density of states at criticalityi 13 : 14 Analogous to 
the conventional quantum Hall effect where the Hall con- 
ductance is quantized, the Hall conductance of the spin 
current of such systems is quantized. This effect is, there- 
fore, called spin quantum Hall effect, because spin, rather 
than charge, is conserved in such a superconductor. 

Studies of the critical behavior of class D with both 
time-reversal and spin-rotation symmetries broken have 
a longer history, starting from the two-dimensional 
random-bond Ising model (RBIM), 17 i 18 i 19 i 20 . 21 i 22 i 23 i 24 
The RBIM can be mapped onto the Cho-Fisher net- 
work model, 19 which resembles the original Chalker- 
Coddington network models for the integer quantum 
Hall plateau transition, but with a distinct symmetry 
The class D models may have an even richer phase di- 
agram; in addition to the spin (or thermal, if the spin- 
rotation symmetry is completely broken) quantum Hall 
phase mentioned above and an insulating phase (which 
is always possible), they may also support a metallic 
phased 

The possibility of spin (or thermal) quantum Hall 
states in unconventional superconductors allows one to 
draw a close analogy between the quantum Hall effect and 
superconductivity, as well as to transfer theoretical or nu- 
merical methods developed in the study of one system to 
another. One of the well-developed numerical methods 
in the study of the quantum Hall effect is the calcula- 
tion of Chern numbers of either single- or many-electron 
states, which allows one to distinguish between current 
carrying and insulating states unambiguously, even in a 
finite-size system. The Chern number method has been 
very successful in the studies of quantum Hall transitions, 
for both integra l 27 i 28 i 29 i 30 and fractional effects j2i as well 
as in bilayer systems^ and also in other contextsiSi 3 ^^ 
For example, by finite-size scaling, the localization length 
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exponent v » 2.3 has been obtainedfS^SLiffli consistent 
with other estimates^ More recent application of this 
method to fractional quantum Hall states (where one in- 
evitably has to deal with interacting electrons) allows 
one to determine the transport gap numerically, 31 which 
is not possible from other known numerical methods. 

In this paper, we report results of numerical studies 
on a lattice model of disordered p-wave superconductors 
with p x +ipy pairing, which conserves the z-component of 
electron spin. As we will show later, this is an example of 
the class D model, and in certain sense the simplest model 
that supports a spin quantum Hall phase. We study the 
localization properties of the quasiparticle states by cal- 
culating the Chern and Thouless numbers of the individ- 
ual states, in ways similar to the corresponding studies 
in the quantum Hall context mentioned above. Physi- 
cally, the Chern and Thouless numbers correspond to the 
Hall and longitudinal spin conductivities in the present 
context, respectively, which are also related to the Hall 
and longitudinal thermal conductivities. We note that 
while it has been pointed out earlier that quasiparticle 
bands or individual quasiparticle states can be labeled 
by their topological Chern numbers^ the present work 
represents the first attempt to calculate them numerically 
and use them to study the localization properties of the 
quasiparticle states in the context of unconventional su- 
perconductors. 

Our main findings are summarized as the following. 
We find that the p-wave model we study supports an in- 
sulating phase and a spin quantum Hall phase with spin 
Hall conductance one in appropriate unit. For relatively 
weak disorder, there exist two critical energies at which 
current-carrying states exist, carrying a total Chern num- 
ber (or spin Hall conductance in proper unit) +1 and — 1, 
respectively; they are responsible for the spin quantum 
Hall phase. Phase transitions between these two phases 
may be induced either by changing the disorder strength, 
or by applying and sweeping a Zeeman field. The field- 
driven transition is found to have the same critical be- 
havior as the integer quantum Hall transition of non- 
interacting electrons. As disorder strength increases the 
two critical energies both move toward E = 0, and an- 
nihilate at certain critical disorder strength, resulting in 
an insulating phase in which all quasiparticle states are 
localized. No metallic phase is found in our model. 

The remainder of the paper is organized as the fol- 
lowing. In section II we introduce the model Hamilto- 
nian of our numerical study, and discuss its symmetry. 
We also describe the application of Chern and Thouless 
number methods to the present problem in some detail. 
We present our numerical results in section III, including 
results of the finite-size scaling analysis of the numeri- 
cal data. Section IV is reserved for a summary and the 
discussion of our results. 



II. MODEL AND NUMERICAL METHODS 

We consider electrons moving on a two-dimensional 
square lattice with linear size L, in the presence of pairing 
and random potentials, described by the Hamiltonian 

n = - f J2 ( c iTS'T + c \i c n) 

<i,3> 
<i,j> 

+ ^(«i-/i)(ct T c iT +ctjc n ), (1) 

i 

where Ui is the random potential on site i which are in- 
dependent random variables distributed uniformly from 
— W/2 to W/2, and [i is the chemical potential of the 
electrons, t is the nearest neighbor hopping integral and 
we choose t = 1 as the unit of energy from now on. We 
define, for a p-wave superconductor, the pairing poten- 
tial Ajj± ex — ±A and Ajj± e = ±iA, where e x and e v 
are unit vectors along x- and y-axis, respectively. For a 
finite-size system, we introduce the generalized periodic 
boundary condition cj+x.^ = e ldi c^, Cj_|_L 4 j_ = e i Cji 
(i = x,y). We note that while the total spin of the system 
is not conserved for the p-wave pairing, the z-component 
of the spin is conserved due to our choice that pairing 
only occurs between electrons with opposite spin. This 
becomes especially clear if we rewrite the Hamiltonian 
in terms of particle-hole transformed operators for the 
electrons with down spins: 

dii = Cj T , dn = ctp (2) 

so that 

W = (44)(i\ 3d, (3) 

where h is the Hamiltonian without pairing potentials for 
each spin component, and A = (A m „). Clearly the num- 
ber of d particles is conserved, reflecting the conservation 
of the z-component of the total electron spin. Thus the 
corresponding transport properties of the z-component 
spin are well-defined; in the following we simply use the 
word spin to refer to its z-component, and spin conduc- 
tances refer to the ratios between the z-componcnt of the 
spin current and the gradient of the z-component of the 
Zeeman field. 

The p-wave pairing symmetry A mn = —A nm leads to 
a special symmetry of the Hamiltonian: 

o x Ho x = -H T , CT * = f i J ) ' ^ 

which makes the current problem distinct from the usual 
problem of electrons moving in a random potential. In 
fact, it is an example of the symmetry class D in the 
classification of Altland and ZirnbaueriS 
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We note that the model we study here, Eq. QJ, i s n °t of 
the most general form of symmetry class T& for spin- 1/2 
fermions, which, by classification, has no spin-rotation 
symmetry along any direction. Consider a generic Hamil- 
tonian for quasiparticles in a superconductor: 



H 



E 

a/3 



(hapC^Cp + ^S a 0clc f fj + ^S* a/3 Cf3C 



(5) 



where a and f} are indices that label both lattice site and 
spin of the electron, running from 1 to 27V (if TV is the 
number of lattice sites). The Hamiltonian can be solved 
by the Bogoliubov transformation, or, explicitly by the 
diagonalization of the 47V x 47V matrix 



H = 



h d 

-6* -h T 



(6) 



In the generic case of class D (without time-reversal and 
spin-rotation symmetries), the only constraint on H is 



W = H 



where 



1 2 N 

Un 



O x ® 1-2N- 



(7) 



(8) 



This constraint, which comes from both hermiticity and 
Fermi statistics, is the same as Eq. (@J). Note that 
this Hamiltonian is twice as large as the p-wave pairing 
Hamiltonian we study [Eq. ©], although they belong to 
the same symmetry class for the following reasons. In- 
terestingly, the partial spin-rotation symmetry along z- 
axis leads to a decomposition of the 47V x 47V matrix H 
into two homomorphic subblocks. One subblock corre- 
sponds to spin-up particles and spin-down holes, and the 
other to spin-up holes and spin-down particles. With- 
out additional symmetry, each subblock belongs to the 
conventional unitary ensemble^ On the other hand, if 
spin-rotation symmetry in other directions is present (as 
in d-wave or other singlet pairing), the spin- up particles 
and spin-down holes are equivalent, and the coupling be- 
tween them is symmetric; this is, in fact, the case of 
class C. In the present case, however, the coupling be- 
tween spin-up particles and spin-down holes is antisym- 
metric, required by the special p-wave pairing we intro- 
duce. Therefore, the model we study is equivalent to 
pairing between spinless or spin-polarized fermions, also 
described by a Hamiltonian of the form [Eq. J^J], with a 
and P label lattice sites only. It is in this sense that we 
can study well-defined spin transport in a class D model. 

It is also useful for us to consider the presence of a 
uniform Zeeman field: 



H B = MS^E( C mT C mT - C lnl c mO 
m 

= HBB ^(<4 T rf mT + 4nl d ml) + c On S t. (9) 



We note that the Zeeman field plays a role of the Fermi 
energy for the (conserved) d particles. More importantly, 
its presence changes the symmetry property of the sys- 
tems, because Hb does not obey Eq. Q. 

The spin Hall conductance of an individual quasiparti- 
cle eigenstate to) can be calculated by the Kubo formula 



(m) 



ih 
~A 



E 



(m\jS\n)(n\jS\m) - (m\j^\n)(n\j^\m) 

En 



(E n — E m ) 2 

(10) 

where A = L 2 is the area of the system. m), n) are 
quasiparticle eigenstates of the Hamiltonian [Eq. (JIJ] and 
Jx i jy the components of the spin current operator. Fol- 
lowing Thouless and co-workers we can show that 
the spin Hall conductance averaged over boundary con- 
ditions is related to a topological quantum number: 



(oL(m)) 



ijj d6 * d6 *li 

I dm dm \ 



dm 
dfh, 



dm \ 
dT x ) 



(11) 



where C\ (to) is an integer and known as the first Chern 
index. As is widely used in quantum Hall transitions 
and other contexts t 2 7i 28 i 29 i 3n i 31 i 32 i 33 i 34 i 3fl d(m) can be 
used to distinguish current carrying states from localized 
states unambiguously, even in finite-size systems, thus 
providing a powerful method to study the localization 
properties of the quasiparticle states. 

An alternative way to study the localization properties 
of the states is to calculate the Thouless number (also 
known as the Thouless conductance) of the states at a 
given Fermi energy E, defined 



9t{E) 



A75 ~ h 1 



(12) 



where AE is the average energy level spacing at energy 
E, and (|5J5|) is the average energy level shift caused by 
the change of the boundary condition from periodic to 
anti-periodic in one spatial direction. It was argued in 
the context of electron localization that gx (E) is propor- 
tional to the longitudinal conductance of the system^i^ 
in the present context we expect it to provide a measure 
of the longitudinal spin conductance of the superconduc- 
tor. Thouless numbers have also been numerically stud- 
ied for the conventional integer quantum Hall effect, in 
both full^ and projected^ lattice models. 

In this work we carry out numerical calculations to di- 
agonalize the Hamiltonian H to obtain the exact quasi- 
particle eigen wave functions. We calculate their Chern 
and Thouless numbers to study their localization prop- 
erties, and perform finite-size scaling analysis to extract 
critical behavior of the transitions driven by the change 
of the disorder strength W or the Zeeman field. 
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III. NUMERICAL RESULTS 



In the absence of the random potential we can diago- 
nalize the Hamiltonian [Eq. l(T]l] in the momentum space, 
and the energy spectrum is 



Eh — 



(13) 



where eu = — 2t(cos k x + cosk y ) — p is the single-particle 
kinetic energy, and A& = 2iA(sinfc :E +isinfc, / ) the p- wave 
pairing order parameter. From Eq. (|13fl . we expect that 
there is an energy gap between two bands, which should 
be stable against weak disorder, while the gap will be 
closed when the disorder becomes strong enough. 

In this work, we choose A = 0.5 and the chemical 
potential p = 3.0 to avoid the van Hove singularity at 
zero energy in the single electron spectrum. To calculate 
the Chern number of each eigenstate, we evaluate the in- 
tegral in Eq. (|llfl numerically over the boundary phase 
space < 9 X , 8 y < 2ir. We divide the boundary phase 
space into M x M square grids with M = 20-80, de- 
pending on the system size L — 10-40 to achieve desired 
precision. Figure shows the density of states (DOS) 
(per lattice site and spin species) p(E) for a system with 
L = 10 and W = 4.0. For such a relatively weak disor- 
der, the superconducting gap is still visible. Also shown 
is the spin Hall conductance function of quasi- 

particle Fermi energy E = psB, calculated by summing 
up Chern number of states below the Fermi energy. We 
find that a xy jumps from zero up by one unit near the 
(disorder-broadened) lower band edge, and jumps back 
to zero above the gap. Therefore, a plateau in a xy is well 
developed around E — 0, clearly indicating the existence 
of a spin quantum Hall phase. This phase with topolog- 
ical Chern number equal to one is the simplest possible 
spin quantum Hall phase for non-interacting quasiparti- 
cles; it is simpler, for example, than the corresponding 
phase of an d x 2_ y 2 + id xy superconductor, which carries 
a total Chern number two. 

In the following discussion, we focus on cases with dis- 
order strong enough to close the gap, and look for tran- 
sitions from the spin quantum Hall phase to other pos- 
sible phases, driven by either the disorder strength W 
or the quasiparticle Fermi energy. In Fig. [3 we plot the 
total DOS p(E) (which is roughly system size indepen- 
dent) and the density of current carrying states (defined 
as states with non-zero Chern number) p e {E) for sys- 
tems with L — 10-40. We find that p e (E) has a weak 
double-peak structure near E = for large L, whose 
width shrinks as L increases. This behavior is reminis- 
cent of those seen in the numerical study of current carry- 
ing states in the integer quantum Hall effectj^&Sii where 
the current carrying states exist only at discrete critical 
energies in the thermodynamic limit and, thus, the width 
of p e {E) peak(s) shrinks to zero as L increases toward in- 
finity. In the present case the two peaks correspond to 
two such critical energies, carrying a total Chern num- 
ber +1 and —1, respectively, which are responsible for 
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FIG. 1: Density of states p(E) (solid line) and spin Hall 
conductance a^ y (E) (dotted line, in units of h/8n) for L = 
10 and W = 4.0. We average over 500 samples of different 
random potential realizations. 
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FIG. 2: Density of states (DOS) p(E) and density of cur- 
rent carrying states (with nonzero Chern number), p e (E), for 
systems with L = 10-40 and W = 8.0. 



the spin quantum Hall plateau when the Fermi energy is 
between them (so that only the lower critical energy is be- 
low the Fermi energy). According to the scaling theory 
of localization, p e (E) depends on L only through a di- 
mensionless ratio L/t;(E) when the system size becomes 
sufficiently large; the localization length £ diverges in the 
vicinity of a critical energy E c as £,(E) ~ \E — E c \~ v . 
Therefore, the number of current carrying states N e (L) 
behaves as 



N e (L) = 2L 2 / Pe (E)dE ~ L 2 - 1 ^, 



(14) 



from which we can estimate v. Assuming we have a simi- 
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FIG. 3: Percentage of current carrying states N e /N(L) versus 
system size L on a log-log scale for W = 8.0. The solid line 
is a power-law fit of the data. 

lar situation here, we plot N e (L), normalized by the total 
number of states N(L) = 2L 2 , on a log-log scale in Fig. [3] 
Just as in the quantum Hall case, 28 - 29 , we can fit the data 
to a power law (a straight line in the log-log plot) as in 
Eq. 114|) reasonably well, and obtain 

v = 2.6 ±0.2. 

This is close to the corresponding exponent v = 2.3 ±0.1 
for the integer quantum Hall transition. These results 
suggest that just as in the case of the integer quan- 
tum Hall effect, current carrying states exist at discrete 
critical energies in the thermodynamic limit, and the 
spin quantum Hall transition driven by the Zeeman field 
(or equivalently, the quasiparticle Fermi energy) has the 
same critical behavior as the integer quantum Hall transi- 
tion. This is expected on the symmetry ground, because 
in this case the critical energies are away from E = 0, 
and thus can only be reached in the presence of the Zee- 
man field. As discussed earlier, the Zeeman field breaks 
the symmetry of Eq. and reduces the symmetry of the 
present problem to that of electrons moving in a magnetic 
field and a random potential. 

While the Chern numbers measure the ability of in- 
dividual states to carry spin Hall current, we have also 
calculated the Thouless conductance gr(E), which is a 
measure of the longitudinal spin conductance. Unlike the 
Chern number calculation which requires the diagonal- 
ization of the Hamiltonian for many different boundary 
conditions, the Thouless number calculation only needs 
the diagonalization at two different boundary conditions, 
thus allowing us to study larger systems. On the other 
hand, it is known in the numerical study of quantum Hall 
effect that Chern number calculation reaches the scal- 
ing behavior at smaller system sizes. Therefore, these 
two methods are complementary to each other. Figure 01 



FIG. 4: Density of states p(E) and Thouless number gr(E) 
for systems with L = 40-80 and W = 8.0. 

shows p(E) and gr(E) for systems with L = 40-80, and 
with W — 8.0. We find that gr(E) has a similar double- 
peak structure as p e (E) with peaks locating at the same 
energies, and that the peaks become narrower as L in- 
creases. In the following we perform the similar scaling 
analysis based on the zeroth moment of gx(E) as for the 
Chern numbers. Namely, we compute the area A(L) un- 
der gr(E) and expect 

/oo 
g T (E)dE ~ IT 1 /". (15) 
- oo 

One slight complication is that unlike p e (E), gr(E) has 
long tails extending to the edges of p(E), which clearly 
has no connection to the critical behavior near the criti- 
cal energies. To eliminate the influence of these artificial 
tails, we introduce a cutoff energy E cut , and exclude con- 
tributions from \E\ > E cut - Based on the Chern number 
calculation above (Fig. |2J), as well as the gr(E) curves 
themselves, we can safely choose E cui between 3.0 and 
4.0, beyond which we find essentially no current carry- 
ing states for L > 40. In Figure [S] we plot, on a log-log 
scale, the area A(L) normalized by the area under the 
DOS curve between — E cut and E cut : 

rE<. nt 

N cnt = / p{E-L)dE 

for a series of different i^cut , and list the corresponding v 
in Table |I| We find that v has very weak dependence on 
the choice of the cutoff energy and its variation between 
2.54 and 2.79 is consistent with the results obtained from 
the Chern number calculation. 

We also studied other disorder strengths. In the case of 
the integer quantum Hall transition^ it is known that 
as the disorder strength increases, the critical energies 
that carry opposite Chern numbers move close together, 
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FIG. 5: Area A(L) of Thouless number gr{E) normalized 
by number of states N cut counted, versus system size L on 
a log-log scale for different cutoff energy iJ C ut and W = 8.0. 
The lines are power-law fits of the data. 
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0.07 



TABLE I: Critical exponent v for different cutoff energy _E cut 
with W = 8.0. 



merge, and disappear at some critical disorder strength 
W c . In the present case, we expect the same to hap- 
pen and due to the symmetry of the Hamiltonian, the 
critical energies can only merge at E — 0. We present 
the results for W = 9.0 in Fig. [H] In this case we no 
longer see two split critical energies, suggesting that the 
two critical energies that were clearly distinguishable at 
W = 8.0 either (i) have moved too close to be distin- 
guishable at the accessible system sizes, or (ii) have just 
merged. We believe scenario (i) is much more likely than 
(ii) based on the following observations, (a) We find that 
the peak value of griE) is independent of system size 
and takes the same value as that of W = 8.0. (b) We 
have performed the same scaling analysis of gr{E) as 
we did above for W — 8.0 and obtained a similar expo- 
nent v « 2.3 (see inset of Fig. |BJ, which is even closer 
to the known value of the integer quantum Hall transi- 
tion. However, there is another possibility that instead 
of entering the insulating phase (in which all quasipar- 
ticles states are localized) immediately, the system is in 
a metallic phase, after the two critical energies merge so 
that the system is no longer in the spin quantum Hall 
phase. Senthil and Fisher— suggested that in this phase 
both the DOS p(E) and the conductance diverge loga- 
rithmically at the band center. Interestingly, we indeed 
find p(E) to be enhanced at E = 0. We believe, how- 
ever, this is not associated to the metallic phase for the 



following reasons, (i) No such enhancement is seen in the 
Thouless number, which is a measure of the longitudinal 
conductance, (ii) We find p(E) to be essentially system 
size independent between L = 40 and L = 80, even at 
E = 0, while one expects^ p(L) ~ logL in the metallic 
phase, (iii) We find that (see below) the enhancement of 
p(E) at E — is also present at stronger disorder when 
the system is clearly insulating. Thus it appears unlikely 
that the metallic phase is responsible for the single peak 
in g T {E). 

The situation is quite different as W further increases. 
In Fig. we present results for W = 10.0 and see a very 
different behavior. Here the peak value of gr(E) system- 
atically decreases as the system size increases, exhibiting 
a characteristic insulating behavior. Combined with re- 
sults of smaller W, we conclude that in the absence of 
the Zeeman field (or when the quasiparticle Fermi energy 
is at E = 0), the system is driven into the insulating 
phase from the spin quantum Hall phase as the disorder 
strength W increases. The critical strength W c is slightly 
above 9.0 and clearly below 10.0. No evidence has been 
found for the existence of an intermediate metallic phase 
that separates these two phases for our choice of model 
parameters (fi — 3.0, A = 0.5, etc.). 

The critical behavior of the transition driven by in- 
creasing W is expected to be different from the one driven 
by changing the Zeeman field discussed above, due to 
the additional symmetry. In order to study the critical 
property one first needs to determine the critical disor- 
der strength W c accurately, which we are unable to do 
within the accessible system size in our study. It would be 
of significant interest to study this transition with more 
powerful computers and/or other computational meth- 
ods. 

We give the results of W — 15.0 in Fig.[S]as an example 
of strong disorder, where all states are clearly localized. 
Here, the Thouless number drops rapidly as the system 
size increases as expected. Interestingly, the enhance- 
ment of the DOS at E = remains to be quite pro- 
nounced, suggesting that it is not associated with pos- 
sible metallic behavior discussed above. For compari- 
son, we have also calculated the DOS for a d x 2_ y 2 +id xy 
superconductor, by choosing the pairing order parame- 
ter to be Ajj +Ex — —Ajj +ey — A x 2_ y 2, Ajj +ex+ey — 
—Ajj +ex - ey = iA xy . In the d-wave superconductor the 
total spin of the system is conserved, due to the sin- 
glet nature of the pairing. As a consequence it belongs 
to the symmetry class C in the classification by Altland 
and Zirnbauer^ This model has been studied in consid- 
erable detail in Refs. lillilllTlTl We plot the DOS for 
a d x 2_ y 2 + id xy superconductor for different values to W 
in Fig. [§] While the gap vanishes just like the p-w&ve 
case for sufficiently large W, the DOS exhibits a pseu- 
dogap behavior at E = for large W, in the vicinity 
of which the DOS vanishes in an (apparently sublinear) 
power law as predicted This is a good example that 
the change of symmetry profoundly affects the critical 
behavior as well as other properties of the system. 
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FIG. 6: Density of states p(E) and Thouless number gr{E) 
for systems with L — 40-80 and W = 9.0. The inset shows the 
area of Thouless number gr(E) divided by iV cu t for E cn t = 
4.0. 



FIG. 8: Density of states p(E) and Thouless number gr(E) 
for systems with L = 20-40 and W = 15.0. 




FIG. 7: Density of states p(E) and Thouless number gr(E) 
for systems with L = 40-80 and W = 10.0. The inset is a 
blow-up of the Thouless number curves near E = 0, which 
shows that gr(E — 0) decreases with increasing L. 

IV. DISCUSSION AND SUMMARY 

In this paper we have studied the localization proper- 
ties of the quasiparticle states in superconductors with 
spontaneously broken time-reversal symmetry, which 
support spin quantum Hall phases. Our study is based on 
the exact diagonalization of microscopic lattice models 
and the consequent numerical calculation of the Chcrn 
and Thouless numbers of the quasiparticle states. Our 
microscopic study is complementary to previous numer- 
ical work on this subject, which have been based almost 
exclusively on effective network models with appropriate 
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FIG. 9: Density of states of d x 2_ y 2 + id xy superconductor 
with L — 60, A ;c 2_ H 2 = 1, A xy — 0.6. We average over 80 
samples of different random potential realizations. 



symmetries, 12 i 14 i 15 i 16 i 19 i 24 

We have focused mostly on a p-wave pairing model in 
which the time-reversal symmetry is broken by the (com- 
plex) pairing order parameter, while the z-component of 
the total spin is conserved so that the transport proper- 
ties of the z-component of the spin is well defined. We 
find the system supports a spin quantum Hall phase with 
spin Hall conductance one in appropriate unit, and an 
insulating phase. Transitions between these two phases 
may be induced either by changing the disorder strength, 
or by applying and sweeping a Zeeman field. The field- 
driven transition is found to have the same critical be- 
havior as the integer quantum Hall transition of non- 
interacting electrons as expected on symmetry grounds. 
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The disorder-driven transition in the absence of the Zee- 
man field is expected to have different critical properties 
due to additional symmetry of the Hamiltonian. How- 
ever, we have not been able to study the critical behavior 
of this transition. 

The symmetry properties of the p-wave pairing model 
in the absence of the Zeeman field belongs to class D in 
the classification of general fermion pairing models of Al- 
tland and Zirnbauer £ It has been suggested that in addi- 
tion to the quantum Hall and the insulating phases, class 
D models may also support a metallic phased which has 
logarithmically divergent density of states and conduc- 
tance. Such a system can have either a direct transition 
between the quantum Hall and the insulating phases, or a 
metallic phase separating these two phases. In our model 
we find a direct transition between the spin quantum Hall 
and insulating phases, but no definitive evidence for a 
metallic phase. This is not unusual as it is known^ that 
specific microscopic models may or may not support the 
metallic phase. 

For comparison, we have also calculated the density 
of states of a d-wave superconductor with d x 2_ y 2 + id xy 
pairing order parameter, which supports a spin quan- 
tum Hall phase with spin Hall conductance two in the 
same unit. This model has different symmetry proper- 
ties and belongs to class C in the classification of Altland 
and Zirnbauer. We find that the density of states van- 
ishes with sublinear power law near E = 0, in agreement 
with earlier studies. 13,14 - 15 - 16 This is in sharp contrast to 
the p-wave case in which we observe an enhanced den- 
sity of states at E = for sufficiently strong disorder, 
demonstrating the profound effect of symmetries on the 



low-energy properties of the system. While this enhance- 
ment is somewhat reminiscent of the divergent density of 
states of the possible metallic phase, further analysis sug- 
gests this is not the case. The origin of this enhancement 
is currently unclear. 

Finally we note that recently there is interest in the 
spin Hall effect in semiconductors with spin-orbit cou- 
pling, which is driven by an electric field4^& Physically 
this effect is quite different from the spin Hall effect dis- 
cussed here, in the following ways, (i) Our spin Hall effect 
is induced by the gradient of a Zeeman field that couples 
to spin, while the other effect is induced by an electric 
field that couples to charge, (ii) The existence of the spin 
Hall effect in our case relies on the broken time-reversal 
symmetry in the pairing Hamiltonian, while the time- 
reversal symmetry is intact in the Hamiltonians used in 
Refs. EHEa: instead the spin-Hall effect is present due to 
the presence of spin-orbit coupling. In Ref. l46l a univer- 
sal spin Hall conductance was found in a clean system; 
it is not clear at present if this value has a topological 
origin as in our case, and how stable this result is in the 
presence of disorder. 
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